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Abstract 


The  paper  shows  that  the  Re. ssner-Mindlln  plate  model  for  clamped-ln 
boundary  condition  does  not  capture  the  boundary  layer  behavior  for  the 
bending  moments.  This  boundary  layer  Is  present  In  the  3  dimensional  formula¬ 
tion.  In  contrast  the  (1,1,2)  model  shows  this  boundary  layer.  The  strength 
of  the  boundary  layer  for  the  (1,1,2)  model  Is  analyzed. 

Key  words:  plates,  Reissner-Mlndlln  model. 
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1. 


Introduction 


There  are  many  plate  models  which  try  to  approximate  the  solution  of  the 
three  dimensional  plate  problem  by  a  system  of  differential  equations  in  two 
dimensions.  Denoting  the  thickness  of  the  plate  by  d,  then  the  difference 
between  the  solution  of  these  models  and  the  three  dimensional  solution  con¬ 
verges  to  zero  (in  the  scaled  energy  norm)  as  d — >0.  The  two  best  known 
models  are  the  Kirchhoff  and  the  Reissner-Mindlin  model.  The  major  difference 
between  various  models  is  the  boundary  layer  behavior  of  the  solution.  The 
boundary  layer  leads  to  large  differences  in  moments  and  shear  forces  in  the 
boundary  region. 

If  the  boundary  is  smooth  and  the  load  is  also  smooth  then  the  Kirchhoff 
solution  is  smooth  up  the  boundary.  In  contrast,  the  Reissner-Mindlin 
solution  has  a  boundary  layer  behavior  whose  strength  depends  on  the  type  of 
boundary  condition.  In  [1]  [2]  a  rigorous  analysis  of  the  boundary  layer 
behavior  for  the  Reissner-Mindlin  solution  on  a  smooth  domain  is  given.  Among 
other  things,  it  is  shown  that  the  boundary  layer  is  strongest  for  the  soft 
simple  support  and  is  very  weak  for  the  clamped-ln  boundary  condition. 

The  Kirchhoff  model,  In  general,  gives  accurate  results  in  comparison 
with  the  3  dimensional  solution  when  the  boundary  layer  is  weak.  Hence  it 
could  be  expected  that  for  the  clamped  boundary  condition,  the  Kirchhoff  as 
well  as  the  Reissner-Mindlin  model  yield  reliable  and  high  quality  results  up 
to  the  boundary.  However,  this  is  not  the  case.  The  Reissner-Mindlin  model 
completely  "misses'4  in  this  case,  the  boundary  layer  which  is  stronger  than 
any  predicted  by  this  model. 

To  show  this,  let  us  consider  the  problem  of  the  uniformly  loaded  square 
plate  w  *  {x^.Xgl  |x^|  <|,  i  ■  1,2)  with  thickness  d  *  1/100.  We  assume 
that  the  material  of  the  plate  is  homogeneous,  isotropic  and  its  Poisson  ratio 
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v  »  0.3.  In  the  Table  1  we  show  the  moments  ^(Xj,0),  M2  2(xj,0)  for 
0<x^<0,5  computed  from  the  3  dimensional  formulation,  Reissner-Mindlin 
solution  (k  =  0.87)  and  the  solution  of  the  model  (1,1,2)  which  will  be 
addressed  in  the  Section  2. 


Table  1.  The  moments  Mj  ^x^.O)  and  M2 


"i.i'V” 

M  , 

2  ,< 

,(xlt0) 

- -  — 

X1 

3D 

RM 

(1,1,2) 

3D 

RM 

1,1,2) 

Ej  l2“/» 

0 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

0 

-0.0229 

0 

.2000 

-0.0157 

-0.0157 

-0.0157 

-0.0163 

-0.0163 

0 

-0.0163 

0 

.4000 

+0.0164 

+0.0163 

+0.0164 

+0.0026 

+0.0027 

0 

+0.0026 

0 

.4900 

0.0470 

0.0470 

0.0470 

0.0141 

0.0141 

0 

0.0141 

0 

.4930 

0.0483 

0.0483 

0.0483 

0.0145 

0.0144 

0 

0.0144 

0 

.4990 

0.0509 

0.0509 

0.0509 

0.0168 

0.0152 

9.2 

0.0170 

1.4 

.4993 

0.0510 

0.0510 

0.0510 

0.0176 

0.0153 

12.3 

0.0179 

2.1 

.4999 

0.0512 

0.0512 

0.0512 

0.0207 

0.0153 

25.3 

0.0212 

2.8 

.5000 

0.0513 

0.0513 

0.0513 

0.0220 

0.0153 

30.0 

0.0220 

0 

We  see  that  on  the  central  line  -  0  the  moment  ^x^.O)  is  very 
accurate,  i.e.  the  error  (with  respect  to  the  three  dimensional  solution)  is 
less  than  0.1%  in  all  reported  points  for  both  the  RM  and  (1,1,2)  models. 

In  addition  the  moment  Mj  ^(Xj.O)  is  smooth  up  to  the  boundary.  In 
contrast,  the  moment  M2  is  smooth  up  to  the  boundary  only  for  the  RM 

model,  and  the  accuracy  is  very  poor  (30%  error).  The  3  dimensional  and 
(1,1,2)  models  clearly  show  the  boundary  layer  and  the  error  of  (1,1,2)  model 
is  practically  acceptable. 

To  understand  the  boundary  layer  better  let  us  define  ^(x^)  so  that 
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exp 


(-e<V 


Hg‘^x1,0)-M2|2(x1.0)| 

|M^M2(0. 5,0)  -  M2  2(0.5,0)| 


The  function  0  characterizes  the  boundary  layer  and  we  can  expect  that  for 
IXj-0.51  small,  0  is  nearly  constant.  M^^tXj.O)  was  used  as  the  (smooth) 
base  function.  In  Table  2  we  show  values  of  0  for  the  3  dimensional  and  the 
(1,1,2)  model. 


Table  2.  The  values  of  the  function  0(x^). 


X1 

3D 

(1,1,2) 

0.4990 

15.43 

13.86 

0.4996 

18.17 

13.41 

0.4999 

23.85 

13.32 

0.4993 

24.13 

13.31 

0.4999 

24.78 

13.18 

We  see  that  the  boundary  layer  of  the  three  dimensional  solution  is  still 
stronger  than  that  of  the  model  (1,1,2);  nevertheless,  as  we  have  seen  in 
Table  1,  the  values  of  the  model  (1,1,2)  seems  to  be  practically  acceptable. 

In  section  3  we  will  show  that  for  the  (1,1,2)  model  we  can  expect 

/120 

0  *  /  l-v  =13.1. 

The  model  (1,1,2)  is  very  robust  for  all  types  of  practically  important 
boundary  conditions.  The  RM  model  can  yield  good  or  bad  results  depending  on 
boundary  conditions.  See  also  [3]  for  further  details. 

In  Section  2  we  will  briefly  mention  the  concept  of  hierarchic  plate 
models.  In  Section  3  we  will  analyze  the  basic  boundary  layer  property  of  the 
model  (1,1,2). 

Computations  reported  in  this  paper  have  been  made  by  the  program 
MSC! PROBE. 
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2.  Hierarchic  plate  models. 

{3  d  3 

x  =  (Xj.x^XjleR  I  (x1,x2)ew,  _2<x3<2|cR  ^  the 

plate  u  with  thickness  d.  By  T  we  denote  the  boundary  of  u.  Further  we 
let 

S  =  |xe  R3|  (Xj.Xg)  €  T,  <  x3  <  | 

R+  =  jxe  R3(Xj  ,x2)  6  T,  x3  = 

We  will  consider  the  plate  problem  as  the  3  dimensional  solution  of 
elasticity  with  homogeneous  isotropic  model  with  |  of  normal  load  on  upper  and 
lower  surface  R+. 

As  usual,  denote  the  displacement  vector  by  u  =  (Uj.u^u^,  <r  =  {o-^  j>, 

e  =  {e,  }  the  stress  and  strain  tensor  respectively  and  by  ^(u)  the 

*■  *  J 

strain  energy  expressed  with  the  usual  Hooke’s  matrix  A  (<r  =  Ae). 

The  exact  3  dimensional  solution  of  the  clamped-ln  plate  problem 
Is  the  minimlzer  of 

GA(u)  =  ^(u)  - Q(u) 
over  all  ueH(Q)  *  {u«  (H1  Cn) )3,  u  »  0  on  S>, 
and 

Q(u)  «  J  §  (u3(xrx.|)  +  u3(xrx2,|))dx. 
a 

By  the  Index  A  we  emphasized  the  use  of  the  Hooke’s  matrix  A.  The  solution 
u  under  our  assumption  satisfies  the  symmetry  conditions:  = 

-u1(x1,x2,-x3),  1  ■  1,2,  u3(xlfx2,x3)  =  u3(x1,x2,-x3) . 

Now  let  n  *  (nltn2,n3),  k 0,  1  ■  1,2,3.  By  the  n-hlerarchical 
solution  u,  we  mean  the  minimlzer  of 


4 


GB(u)  *  ^(u)  -Q(u) 


over  the  set 

(2.1) 


X^n'(Q)cJ((0)  of  all  functions  of  the  form 

ni 

Uin)(xl'X2’X3)  "  ^  Uinj^Xl,X2^X3’  1  =  1'2,3- 
J=0 


B  denotes  a  Hooke’s  matrix,  possibly  with  modified  coefficients.  Under 

our  assumption  the  symmetry  of  the  load  yields  u,  =  0  for  J  even  and 

1  *  J 

i  a  1,2  and  u_  .  ■  0  for  J  odd. 

J 

We  have  shown  in  [4]  that  the  solution  of  the  model  (1,1,0)  is  identical 
with  the  Reissner-Mindlin  solution  provided  a  modified  matrix  B  (B*A)  is 
used.  The  solution  of  the  model  (1,1,2)  mentioned  in  the  Table  1  is  computed 
for  B  =  A. 

We  defined  the  solution  for  clamped-in  plate  only.  It  is  obvious  that 
for  general  homogeneous  boundary  conditions  only  the  constraint  of  K(Q)  on 
S  has  to  be  modified. 

If  the  solution  u  Is  independent  of  and  u2  =  0,  then  the  Reissner- 

Mindlin  model  becomes  the  Timoshenko  beam  and  {u^n^>  becomes  a  hierarchy  of 
the  beam  solutions.  From  (2.1),  we  can  compute  the  moments  M  ,,  and  M0  0  in 

M  t  «- 

the  usual  way. 

3)  The  hierarchy  of  the  models. 

In  this  section  we  will  consider  the  hierarchy  of  the  beam  models  for 
plain  strain.  This  is  the  special  case  of  the  plate  when  the  solution  is 
independent  of  x^  and  -  0  (for  example,  an  infinite  strip  plate). 

The  model  (1,1,0)  (with  modified  matrix  B)  leads  to  the  well  known 
Timoshenko  beam  equation  whose  solution  has  no  boundary  layer. 

The  model  (1,1,2)  uses  the  set  of  functions 
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Ui<x1,x3)  -  ut  l  1x^X3 


u3<xl.x3)  '  U3,0(xl’*U3,2(xl)x3 

with  u 2  =  0.  The  Euler  equations  for  the  functions  u^  u^  q,  u^  2*  °f  the 
minlmizer  of  GA  are 

-12  j  *  (2X-M)u^  2  =  0 

(3.1) 

-12  Md'2!--,  ,*<>3  „)-  -  ^2  =  12<! 

4(A+2(i)u3, 2  -  2XuJ  j-lil-u,  „)' 

3  .2  „  ~,2 

~  20^  U3,2  =  341  q 

where 

„  Ev  E 

A  =  (1+v) (l-2v)  ’  M  2(l+w)  ‘ 


The  boundary  conditions  for  the  clamped  beam  are  Uj  1  s  “3  q  s  “3  2  = 
System  3.1  can  be  solved  explicitly.  Assuming  for  simplicity  that 
-  | <  xi < |  and  q  is  symmetric  then  the  solution  is 


f  f 

x,  —  x. 


,  x  fd  1  d  ll 

ul,l<xl’  “  cl[e  -•  j 


v'So  /P5 


a  +  c. 


|W  1  1 

l"2E  ri2XlJ 


v(l+v) 

2E 


[l  + Yojd2  ed  1  e  q(y)dy+  e  d  1 


e  _x,  r 

1  -y 


e  q(y)dy 


i>(l+v) 


(‘’ToKJ 


X1  X1 

,  , .  12(l+i>2)  f  (x!-y)2  .  ,  . 

q(y)dy  +  - gg -  — ^ —  q(y)dy 
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(3.2) 


f  f 

-x.  — x 


w4il-il)  srd2,c2ifivc: 


.  (l+v)(v+10)(5-v)v'T:i>  j2 
+  — — - -  cl 

40v'30  E 


rt'  ~AX' 

ed  e  g(y)dy  -  e  d 
0 


r 


1  <? 

e  g(y)dy 


1 +v>  2  2 

+  TO?  {v +5v~20)d* 


'  ,  ,  ..  12C1-IT) 

(Xj-yJgtyJdy  +  - g - 


kl 


(xr-yT 


-g  ( y )  dy 


f 

— x. 


,  ,  „  fd  1  d  ll  „  v(l+v)  1 

U3,2  X1  '  Cl[®  6  J  C2  4E  T2 


where 


3( 1-v)  VY^v 
L  2/30 


♦  6 


w(l+w) 


(v+10)d 

r  f 

5Xl 

e 

-X,  f  f 

1  'Sy  ';Xl 

e  g(y)dy  -  e  d 

1  5y  1 

e  g(y)dy 

,  « 

. 

J 

0  0 


(x1-y)g(y)dy 


0 


f 


and  C^,  C^,  are  constants  which  are  determined  from  the  boundary  condition 
at  Xj  =  |  (because  of  symmetry).  Formulae  (3.2)  show  that  the  boundary 
layer  for  the  model  (1,1,2)  for  the  beam  is  of  order  v  provided  that  the 
constant  Cj*0.  In  the  case  of  the  clamped  boundary  for  w  >  0  In  fact 
Cj^O.  The  moments  Mj  jfXj)  and  M2  2^X1^  are  comPuted  from  u^  j^Xj). 

In  the  similar  way  also  higher  models  could  be  investigated. 

The  behavior  of  the  moment  2^xi’®^  mentioned  in  the  Section  1  is 
very  close  the  the  beam  behavior.  In  fact  we  see  In  the  Table  2  very  good 
agreement  with  the  strength  of  the  boundary  layer  described  here. 

Finally  we  show  another  simple  numerical  example.  Consider  the  plane 
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strain  elasticity  problem  on  the  domain  0  »  {IXjl  <  .5,  |Xg|  £  0.005}  vlth 
clamped-ln  boundary  conditions  for  |x^|  =  0.5. 

7 

Consider  the  case  when  E  =  10  ,  v  =  0.3  and  a  uniform  load  is  imposed  on 
the  upper  and  lower  side.  Because  of  symmetry,  only  one  half  (x^  >0)  of 
the  domain  will  be  considered.  The  solution  of  the  model  (1,1,2)  Is  solved 
by  the  p-version  of  the  finite  element  using  only  2  elements  (0,0.492), 

(0.492,  0.5)  of  degree  p.  For  the  analysis  of  the  p-version  for  the  beam 
problem  with  relation  to  locking  effects,  we  refer  the  reader  to  [5]. 

In  Table  3,  we  show  the  strain  energy  £  of  the  finite  element  solution 
as  a  function  of  p. 

Table  3.  Energy  6  of  the  finite  element  solution  as  function  of  p. 


p 

g.103 

2 

0.313411 

3 

0 . 376942 

4 

0 . 378874 

5 

0 . 379102 

6 

0.379162 

7 

0.379176 

8 

0.379178 

Table  4  shows  the  moment  Mg  gCx^)  ^or  the  beam  understood  as  special  case  of 
the  plate)  as  function  x^  computed  by  elements  of  degree  8. 

Table  4.  The  moment  Mg  g(Xj)  ^or  the  (1***2)  beam  model 
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X1 

M2,2fxl)1°2 

0 

-1.2496 

0.2 

-0.4966 

0.4 

+1.1503 

0.45 

1 . 7878 

0.495 

2.4225 

0.4977 

2.5216 

0.4995 

3.0536 

0.49995 

3.5062 

0.5 

3.5723 

Approxlmat ing 


-13.lf0S~Xll 

i0\ ,2<*i>  ■»'  **Be  I  d  >  =  H2  2(x, 


)10^  with  A  = 


2.4673641,  B  =  1.1049693  we  get  the  values  reported  in  Table  5.  We  also 
write  the  relative  error  in  percent  when  values  from  Table  4  are  used. 


Table  5.  The  approximate  moment  2^xi^  its  error  e 


X1 

fi2,2  102 

e  */. 

0.495 

2.4689 

1.9 

0.4977 

2.5216 

0 

0.4995 

3.0413 

0.4 

0.4995 

3.5022 

0.1 

0.5000 

3.5723 

0 

We  see  that  the  values  in  the  table  have  the  predicted  boundary  layer 
behavior.  We  have  very  different  behavior  of  Mj  jCXj)  and  M ^  2^xl^'  The 
equilibrium  condition  directly  shows  that  the  moment  Mj  ^  cannot  have  a 
boundary  layer.  On  the  other  hand  the  boundary  layer  of  displacement  vectors 
has  to  lead  to  a  boundary  layer  of  M,  0- 
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Hie  Laboratory  for  Numerical  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Maryland,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals: 

•  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  m  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  die  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  indudes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  L  Bahuaka,  Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Sdence  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742. 


